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A method is presented, which employs the density matrix renormalization group technique in order 
to construct exact ground state exchange correlation functionals for models of correlated electron 
systems coupled to external reservoirs. The technique is applied to the M-site resonant level model. 
We calculate its exact Kubo-conductance, which is available within DMRG, and compare to the 
single-particle conductance obtained from Kohn-Sham energies and orbitals of the exact ground 
state density functional theory (DFT). It is found that the position of transport resonances is 
reproduced essentially exactly, while deviations in the level broadening can be less than 1% and do 
not exceed 10%. Our findings lend strong support to a recently held point of view, namely that 
approximations in the ground state functionals used in DFT based transport calculations can lead 
to drastic errors in transport calculations while exchange correlation contributions to the induced 
effective potential tend to be less significant. 
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Running an electrical current through individual 
molecules and being able to control the current flow by 
molecular design is one of the intriguing aspects of Molec- 
ular Electronics. Naturally, the interesting behavior of 
such "devices" is frequently related to specific properties 
of the molecule at hand. Therefore, in order to prop- 
erly understand experimental findings, notably the cur- 
rent voltage characteristics, ab initio calculations play a 
key role in this field. 

The challenge in such ab initio calculations for trans- 
port properties is to accurately describe interaction ef- 
fects on the junction together with wavefunction hy- 
bridization with the contacts. To deal with the latter, 
it is mandatory to treat a sizable part of the contact on 
the same footing as the molecule. To achieve sufficient 
system sizes, the current "standard method" is based on 
the density functional theory (DFT).[l|,|3,[l] 

A justification of the principal approach, namely the 
use of the single particle Keldysh formalism in conjunc- 
tion with the Kohn-Sham (KS) scheme of time dependent 
(TD) DFT, has been proposed in Ref. on a heuristic 
level and by Stefanucci and Almbladth has been based 
on a detailed analysis using the Keldysh-technique.Q. 

In practical calculations approximate exchange corre- 
lation (XC) functionals have to be employed. Conse- 
quences of such approximations have been investigated 
in Ref. where a relation between the linear current 
j, the KS-Kubo-conductivity and the TDDFT functional 
Exclri^) has been derived: 
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(Etot(rw) denotes the sum of external and Hartree field, 
no(r) the ground state electron density.) In one dimen- 
sion and the dc-limit, j(r, w 0)=/, so Eq. ([T]) reduces 
to / = (^tot + ^c) introducing the KS-conductance, 
gs, the physical voltage Vtot and its XC-shift Vyic [I]- 



In most current implementations of the standard 
method, the linear current response derives from cal- 
culating in a local density approximation. More- 
over, the XC-terms, Vxc or Exc, are universally ignored 
even though Exc(t^) is generally important for the cal- 
culation of excitation energies [6|. Such approximations 
must be expected to have severe consequences. Most 
notably, the missing derivative discontinuity in avail- 
able XC-functionals can artificially boost the DFT level 
broadening Fdft observed in the variation of Qs with the 
gate voltage, V^atc , from its true value F up to the inter- 
action energy [/ 3> F and thus completely impair an ac- 
curate prediction of transport coefficients [sj . This effect 
is a major suspect to cause much of the strong deviations 
that are frequently observed between theoretical and ex- 
perimentally determined transport coefficients. (sl. [ol. [lot 

Further analysis of Eq. ([1]) has been impaired by the 
fact, that for correlated molecular systems, which have 
been coupled to external reservoirs, exact XC function- 
als are not available. To overcome this fundamental bar- 
rier, we put forward a general idea in this work, namely 
to marry DFT with the density matrix renormalization 
group (DMRG) technique. Our approach allows to con- 
struct exact ground state (GS) functionals for a broad 
class of generic model systems of correlated fermions in- 
cluding external reservoirs. Here, we apply this method 
to obtain the exact KS conductance gs ■ By comparing to 
the exact physical conductance, g = I /Vtot, it is demon- 
strated that (7s gives a highly accurate estimate for po- 
sition and broadening of transport resonances implying 
that dynamical corrections Vxci including the viscous 



part [ill, are small. We show this specifically for the in- 



teracting resonant level model (IRLM) and the resonant 
chain model under conditions of the Coulomb blockade. 
Arguments are given, why we expect our results to be 
more generally valid. 

The IRLM and its extension including reservoirs, that 
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FIG. 1: (Color online) Schematic representation of the calcu- 
lational setup. Filled circles indicate the extended molecule, 
eM. Gate voltage, V^atc, is applied to molecule (dot), M, 
only. Interacting sites are coupled by fat line links. Link labels 
indicate hopping amplitudes in the leads t, on the molecule 
tdot and near the contacts t' and t^. (t=l, tdot=0.5 in all cal- 
culations presented.) TZljTZr denotes the reservoirs, here in 
fc-space representation. Upper panel: single site resonance 
level model. Lower panel: five site model. 



we employ in this work, is a suitable testbed to study 
strong correlation effects on transport. The model 
Hamiltonian reads Ti. = TleM+'^u +^tz+'^t , with 
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where c| (c^) and (c^.) are spinless fermionic creation 
and annihilation operators at site £ (at momentum fc), 
■hi = c\ci. Furthermore, 'Hn=J2keiZL,nB. '^kc\c^. and 



( E ^fc^ci + E ^kcicM,) + h.c. (4) 
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Indices denote Hilbert spaces of the molecule proper 
(nanostructure/quantum dot), Ai, of the extended 
molecule, eM and of the left and right electrode reser- 
voirs, TZl^TZr. TCt denotes the tunneling Hamiltonian 
describing the contact between reservoirs and the ex- 
tended molecule; the interaction on the molecule is mod- 
eled by Hu- {U=2t in all calculations presented.) The 
indices 1 and Me denote the first and last site in eM. 
The general setup is displayed in Fig. [TJ For the single 
level model the interaction U extends on the link between 
the resonant level and the neighboring lead sites. In all 
calculations we have considered the case of a half filled 
band and zero temperature, £'Fcrmi=0,T = 0. 

We calculate the GS of our model Eg. ([2]) by means 
of a DMRG calculation. DMRG[i3, is a method 
that searches for an optimized subspace of the complete 
Hilbert space in which selected many body states can be 
described accurately. Notice, that our setup include the 
leads in momentum space representation, Ti-jz, which is 
non-standard, but crucial for later transport calculations. 
Fig. [TJ for details see. ij] In this way we obtain the lo- 
cal electron density on eAi together with the occupation 



number of lead levels. Typically, we use more than 1300 
states per block and ten finite lattice sweeps. 

Next, we briefly explain how an exact DFT may be 
constructed generalizing earlier ideas by Gunnarson and 
Schonhammer. [is} We define a KS-Hamiltonian Hs = 
Ti-Q + Vhxc, with Ho—HeM+'^n+HT (free fermions) 



and an XC potential Vhxc=X], 



also including the 



Hartree term. Here, rij denotes the particle density. The 
sum is over the lattice sites i of eM and the states k of 72.. 
According to theorems by Kohn and Hohenberg[3| and 
Kohn and Sham[17[ a unique set of coefficients Vj speci- 
fying the XC kernel exists, such that the KS-particle den- 
sity coincides with the exact density of the many body 
GS. In practice, we find vj solving a standard optimiza- 
tion problem, which yields a final relative density mis- 
match of less than 10"-'^° per site. 

Within DMRG the linear conductance can be obtained 
from evaluating the Kubo-formula 0,11 
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where is the many body ground state, ry is the broad- 
ening parameter. J„ is the current density operator at 
the bond between site n and n — 1 and Eq denotes the 
GS energy. Due to particle number conservation, the dc- 
conductance is independent of ni,rt2. For the details of 
the procedure see Ref. 
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The conductance (7 - as in fact any dynamical corre- 
lator at T=0 - can be calculated evaluating proper GS 
matrix elements of certain known many body operators, 
see e. g. Eq. Therefore, the general principles of 

DFT apply and functionals exists, parameterized by 77 
(or frequency w, times t, t' etc.), which allow to calculate 
such correlators from the GS density no alone. The func- 
tional, which would yield the exact conductance ^[no] is 
not known. An approximative expression for g coincid- 
ing with the exact KS conductance Qs is readily obtained 
employing Eq. ([5]) using the KS ground state with KS 
single particle energies e^, the corresponding orbitals for 
evaluating the matrix elements Jop and Ti^TLs 
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/(e) denotes the Fermi-Dirac occupation numbers. 

We have calculated conductances for molecules with 
only one site (single interacting level) and with five sites. 
We begin with the single level model with numerical 
parameters: t' ~ 0.1. The model enjoys a particle 
hole symmetry, so that the single transport resonance is 
pinned to the band center, E—Q. Fig. [2]shows the exact, 
g{Vga.tc)^ and the Kohn-Sham conductance gs- Compari- 
son to the non-interacting limit (U—Q)) exhibits a strong 
(280%) interaction driven enhancement of the resonance 
width, r — 0.116, compared to the non-interacting case 
r = = 0.04. Exact DFT, gs, is able to reproduce this 
renormalization effect with accuracy better than 10%, 

rDFT=o.io6. 
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FIG. 2: (Color online) Linear conductance over the gate volt- 
age for the IRLM model with t' = 0.1. Comparison of the 
conductance calculated within DMRG for the full problem 
(+) and for the corresponding effective DFT Hamiltonian 
(x). The line through the DMRG data is a guide to the eyes 
and the Lorentzian is the non-interacting result as reference. 
DMRG half width: V = 0.116 ± 0.001; DFT: 0.106 ± 0.002; 
noninteracting system, f/ = 0: 0.04. 



FIG. 3: (Color online) Comparison of the exact conductance 
(-I-) and the ground state DFT approximation (x, dashed line 
as a guide) for a five site system (t'=0.2). For comparison the 
noninteracting f/=0) is shown as well (long dashed line) . Solid 
line indicates the particle number A'^^iCl^atc) the molecule. 
The resonance of g are sitting Vgato = 0, 1.854, and 2.779 
with a resonance width of F = 0.026, 0.015, and 0.0033. 



We now turn to the five site case, which afFords four 
single levels, that are not pinned to zero energy. Since 
now resonances experience an occupation dependent and 
interaction driven shift ( "Coulomb blockade" ) with a cor- 
responding change of the resonance width, this model can 
serve to investigate the DFT handling of such renormal- 
ization phenomenona. 

Fig. [3] displays the ^(Vgatc) and gs(V^ato) conductance. 
Since particle-hole symmetry implies invariance under 
14;ate — ^^atoj Only the positive branch is shown. The 
first resonance at nonzero energy signalizes the transi- 
tion, where the electron number N M^^a.ta) of the molec- 
ular dot changes between two and one, see Fig. [3l This 
happens at VJjatc ~ U—A, where A denotes the single 
particle level spacing. This expectation is roughly con- 
sistent with the numerical value 1.8 obtained from Fig.[3l 
for details see [l8|. 

Fig. [3] clearly shows, that the DFT calculation per- 
fectly well captures the position of the transport reso- 
nances. In addition, the broadening of the resonance 
peaks is described reasonably well. Similar to the single 
level case, for the center peak width 10% deviations have 
to be accounted for. By contrast, a logarithmic plotting 
is required in order to make the relative deviations visi- 
ble for the broadening of the shifted peaks. Remarkably, 
near resonances the variation of the conductivity is de- 
scribed over more than three orders of magnitude with 
deviations of a few percent or less. 

In Fig. [4] we show the evolution of the local on site 
potential u^^'^ with increasing gate voltage. The overall 
behavior is complicated, and a detailed discussion has 
to be relegated to Ref. [l3|. Here, we can only briefly 
comment on two crucial aspects. First, w^^^ partially 
compensates Vgatc for repulsive voltages inbetween two 
resonances, keeping Nm integer, c. f. Fig. [3l Second, 
at the resonance, V^atc~l-855, the center peak of Fig. [4] 



rapidly decays. Thus, the double well structure in the full 
effective potential v^'^'~^+Vgatc, that was appropriate for 
two repulsive particles Nm='2, transmutes into a single 
well hosting the lone particle, Nm=^- 

The predictive power of conductance calculations with 
ground state DFT may seem surprising at first sight, be- 
cause Vxc is neglected albeit it is very well known that 
the bare KS-response yields incorrect excitation energies, 
which are shifted to proper values in TDDFT only by in- 
cluding dynamical correlations. [20j| A closely related fact: 
the bare KS-spectral function, As(uj,Vgi,to = 0) exhibits 
excitation peaks at frequencies uj of the order of A, and 
not U. The point, that we wish to make here is, that 
for correlated electron systems the dependency of on 
its arguments uj and V^ato is quite different. The linear 
transport probes only in the vicinity of zero frequency. 
The evolution of As{uj = 0, V^atc) with gate voltage is 
closely tied to the particle number and therefore can be 
physically meaningful and give quantitative results even 
if dynamical corrections are ignored. 

We elaborate on this issue and give two reasons, why 
in fact the good performance of is not entirely unex- 
pected. (1) We are concerned with isolated resonances 
r ^ A, J7, which is a situation typical of the relatively 
small molecules that one deals with in the field of Molec- 
ular Electronics. These resonances occur precisely at the 
degeneracy point, where the N and A^+l particle states 
of the molecular dot coincide in energy, so that the parti- 
cle number is half integer N+1/2. Since the exact DFT 
monitors the true particle number on the dot, the de- 
generacy point of the KS-occupation and hence the KS- 
transport resonance coincide with the true value. (2) It is 
much less obvious, why also the resonance width F should 
always be given very accurately, and in fact there is no 
reason to believe that this is the case. However, under 
fairly general assumptions, one may argue that the width 
31] of the transition in NMiVgntc) gives a very good es- 
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to be addressed in future work, to what extend DFT 
can capture also those phenomena - be it in the density 
response or in transport signatures - which hve beyond 
the regime of attraction of the Fermi hquid fixed point. 

Useful discussions with K. Burke, F. Furche, G. Schnei- 
der and P. Wolfle are greatfuUy acknowledged. This work 
was supported by the Center of Functional Nanostruc- 
tures at Karlsruhe University. 



X 



FIG. 4: (Color online) Potential vf^'^ on sites £ of eM cor- 
responding to conductance data shown in previous Fig. |3] 
at Vgate=0.2, 0.005, 1.5, 1.85, 1.86, 2.5, 2.8 (labeUng center site, 
x—6, from top to bottom). 



timate for the width F of the transport resonance. In- 
deed, consider a molecular dot invariant under exchange 
of left and right reservoirs (symmetric coupling). Then 
the single particle lifetime proper, of molecular 

excitations also appears as a transport rate F w QE, 
since the escape rates into the left and right leads, F^ 
and Ffl, simply coincide with 3S/2. Therefore F~^ sets 
the (only) time scale for relaxation processes and hence 
it should also describe the width of transport resonances. 
Notice, that A0\4(V^atc) cannot easily distinguish two sit- 
uations, where F^ and Tn are vastly different with the 
sum, F=Fl-|-F/j, being kept fixed. The excellent agree- 
ment found in this work for the symmetric case may no 
longer pertain into the strongly asymmetric limit, since 
gs (X FiFfl. 

In summary, we have presented a method for perform- 
ing exact DFT calculations for model systems based on 
the density matrix renormalization group (DMRG). The 
approach has been used in order to calculate ground 
state Kohn-Sham conductances for the interacting res- 
onant level model (IRLM), which can be compared to 
exact resuhs obtained with DMRG. We find that DFT 
calculations can describe positions and broadenings of 
transport resonances with a very good accuracy. In fact, 
the resonance position will be given very precisely for any 
correlated electron system, as long as it is connected to a 
single resonant free fermion level by adiabatically switch- 
ing on the interaction (Fermi liquid regime). A further 
implication suggested by our result is that dynamical cor- 
rections should be small as long as vertex corrections can 
be ignored. So the most pressing limitations in practical 
DFT conductance calculations appears to be the missing 
deriviative discontinuity. This poses a notoriously diffi- 
cult problem which, however, in principle has been well 
understood. 21 1 

Finally, we mention that our method DFT®DMRG is 
a very general approach, and not restricted to the IRLM 
employed in this work. It will an intriguing question 
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